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The master equation for a probability density function (pdf) driven by Levy noise, if conditioned 
to conform with the principle of detailed balance, admits a transformation to a contractive strongly 
continuous semigroup dynamics. Given a priori a functional form of the semigroup potential, we 
address the ground-state reconstruction problem for generic Levy-stable semigroups, for all values 
of the stability index /i £ (0, 2). That is known to resolve an invariant pdf for confined Levy flights 
(e.g. the former jump-type process). Jeopardies of the procedure are discussed, with a focus on: (i) 
when an invariant pdf actually is an asymptotic one, (ii) subtleties of the pdf ^-dependence in the 
vicinity and sharply at the boundaries and 2 of the stability interval, where jump-type scenarios 
cease to be valid. 

I. CONCEPTUAL BACKGROUND 

As the Levy noise model we consider any member of a subclass of uni-variate stable probability distributions 
determined by a characteristic exponent —F(p) = — \p\^ of < exp(ipA) >, with < /i < 2. The induced free jump- 
type dynamics, < cxp(ipX t ) >= cxp[— tF(p)], where t > 0, is conventionally interpreted in terms of Levy nights 
and quantified by means of a pseudo-differential (fractional) equation for a corresponding time-dependent probability 
density function (pdf) 

d t p = -\A\^ 2 p = ( [w^x\y)p(y) - w^y\x)p(x)]dy . (1) 



The jump rate w fJl (x\y) cx l/\x — y\ 1+ ^ is a symmetric function, w^(x\y) — Wft(y\x). 

The "free" fractional Fokkcr-Plack equation ([I} has no stationary solutions. On the other hand, asymptotic invariant 
pdfs for confined Levy flights are known to arise in the standard Langevin modeling of an external forces impact on 
the Levy-stable noise, [l| . A disadvantage of that approach is a non-existence of Boltzmann-type (thermal) equilibria, 
0-0 ■ We note in passing that the reference stable laws generically have no moments of order higher than one (and 
may as well have none at all). To the contrary, pdfs for confined Levy flights may have an arbitrary, not necessarily 
finite number of moments. 

Interestingly, if we enforce [1,0 the principle of detailed balance to hold true, by a suitable modification of transition 
rates w(x\y) wu(x\y) (with U playing the role of an external microscopic potential), asymptotic Boltzmann-type 
equilibria of the form p{x) ~ exp[— U(x)] in principle become admissible. The price paid is that the standard Langevin 
modeling of confined Levy flights becomes inadequate. 

We generalize the master equation ([T]) to encompass non-symmetric jump rates as follows, w^(x\y) — > w^(x\y) ^ 

W V(x\y) = wMv) ex P ( U{V) - U{X) ) , (2) 



where U(x) is a continuous function on R. With w^(x\y) replacing w^(x\y), Eq. (JTJ takes the form 

a t(0 = -[exp(-[//2)]|Ar/ 2 [cxp((7/2) (0 ]+pcxp((7/2)|Ar/ 2 cxp(-C//2). (3) 

For a suitable (to secure normalization) choice of U(x), p e q(x) oc exp[— U(x)} is a stationary solution of Eq. ([3]). The 
detailed balance principle necessarily follows: 

wu{x\y)p eq {y) = wu(y\x)p eq (x) . (4) 

The master equation (J3|) cannot be derived within the multiplicative or additive Langevin modeling. 

As mentioned before, Eq. ([3]) at least on formal grounds, admits a transformation to a strongly continuous semigroup 
dynamics. That is akin to a mapping of the standard Fokker-Planck equation into the generalized diffusion (semigroup) 
equation which is widely exploited in the context of diffusion-type processes (there e.g. the Fokker-Planck operator is 
mapped into a symmetric operator, which furthermore needs to be extended to a self-adjoint one). In the Brownian 
case, the Langevin and semigroup derivations, while interpreted in terms of p(x,t), refer to the same diffusion-type 
process. To the contrary, this equivalence does not persist in case of confined Levy flights. 
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The passage from Eq. ([3]) for confined Levy flights to the semigroup dynamics is accomplished by means of a 
redefinition 

p(x,t) = pl /2 (x)V(x,t) (5) 

where p*(x) = p eq {x) = Z^ 1 exp[— U{x)} is an asymptotic invariant pdf of the jump-type process, while the dynamics 
of a real positive-definite function \P(x, i) follows the semigroup pattern [exp(— Ht) ^](x, 0) = ^(x,t) for t > 0, c.f. 
0-0]. Here, we have introduced the Levy-Schrodinger Hamiltonian operator with an external potential 

H^ = \A\^ 2 +V(x). (6) 

Suitable properties of V need to be assumed, so that (i) — is a legitimate (self-adjoint) generator of a (strongly 

1/2 

continuous, contractive) semigroup exp(— tH^), (ii) asymptotically, as t — !> oo we get ty(x, t) — > ^*(x) ~ pj (x) and 
'J'* (2) is a unique ground state of H. Said otherwise, once we have a priori selected an invariant probability density 
p eq (x) = p*(x) oc exp[— C7(x)], ([3]) to justify its interpretation as an asymptotic pdf of a well defined jump-type process 

1 /2 

we must have guaranteed an existence of an associated contractive semigroup dynamics for which ty*(x) ~ pj (x). 
Looking for stationary solutions of the semigroup equation dt^f = —H^, we realize that if a square root of a 

1 /2 

positive invariant pdf p it {x) is asymptotically to come out via the dynamics — > pj , then the resulting fractional 

1/2 

Sturm-Liouville equation H^pJ = stands for a compatibility condition upon the functional form of V(x): 

IAKV /2 

p* 

We note that the inferred semigroup dynamics provides a solution for the Levy stable targeting problem, with a 

1 /2 

predefined invariant pdf. It is an identification p[xA) = pj (x) 1 i(x,t) that does the job. We have discussed this issue 



in some detail in our previous publications, [H, 

Inversely, if we choose a priori a concrete potential function V(x), then an ultimate functional form of an invariant 

1 /2 

pdf p*(x) (actually pj {x)), if in existence, needs to come out from the above compatibility condition. However, a 
solvability of ([7]) with respect to pj merely associates an invariant pdf with a pre-defined semigroup potential V. 
In regard to the asymptotic regime ^(x,t) — > p X J 2 we need much more. Namely, p 1 J 2 {x) must belong to the domain 
of a self-adjoint operator H and a proper identification of such domain may sometimes become tricky. 

The latter problem, and jeopardies involved, we address in the present paper for a predefined function V(x), while 
admitting all stability index values < p < 2 in the compatibility condition ([7]). 

Remark 1: For clarity of presentation and self-explanatory features of discussion we add few comments about the 
stochastic process in question. Let H be a self-adjoint operator in a suitable Hilbcrt space domain. Additionally let 
V = V(x) be a bounded from below continuous function. Then, the integral kernel k(y, s, x, t) = {cxp[— (t— s)H]}(y, x), 
s < t, of the semigroup operator exp(— tH) is positive and jointly continuous in all variables. The semigroup dy- 
namics reads: ty(x,t) = J ^(y, s) k(y, s, x, t) dy so that for all < s < t we can reproduce the dynamical pattern 
of behavior, actually set by Eq. ([3]), but now in terms of Markovian transition probability densities p(x, s,y,t): 

p(x,t) = p*(x)ty(x,t) = J p(y,s,x,t)p(y,s)dy, where p(y, s, x, t) = k(y,s,x,t) pj {x)/pj (y). An asymptotic be- 
1/2 

havior of ty(x,t) — > pj (x) implies p(x,t) — > p*(x) as t — > 00. 

Remark 2: The spectral theory of fractional operators of the form © has received a broad coverage in the 
mathematical [8l4l2l| and mathematical physics literature [ill 13 • Various rigorous estimates pertaining to the decay 
of the eigenfunctions at spatial infinities, quantify the number of moments of the associated pdfs for different classes 
of potential functions V(x). As well, fractional versions of the Feynman-Kac formula for an integral kernel of the 
semigroup operator have an ample coverage therein. 



II. /i-FAMILY OF SEMIGROUP POTENTIALS FOR A PREDEFINED INVARIANT PDF AND THE 

p G (0, 2) BOUNDARY ISSUE 

For a pseudo-differential operator |A| M / 2 , the action on a function from its domain is greatly simplified in the 
physics-oriented research. Normally, with v^(dx) standing for the Levy measure, we have: 

(|Ar /2 /)(z) = - / [f(x + y)- f(x) - u^dy) (8) 

Jr 1 t y 
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By turning over to the Cauchy principal value of the involved integral, one actually considers 

(\A\^f)(x) = - J [f(x +y)- f(x)]^(dy) (9) 

which upon changing an integration variable y — > z = x + y, may be reproduced in form, see e.g. Q 

_ (|Ar/2 m = I> + l)sin (7 r M /2) I dz {l0) 

7T J \Z X\ 

Let us investigate the properties of the — | A|^/ 2 /(a;) by turning over to the Fourier image of f(x). We employ a 
redefinition of Eq. (fTU| 



which yields 



|Ar / Vw = £<i±M r /(t)e ..^ r tr^. (12) 

IT'S/ Z7T J — oo J — co 12/1 



The integral over can be calculated as follows 



(e-*-l)cfo _„ r (ooBfcy-1)^ _ Wr( _ M)cog Z[M. (13) 



|y|l+M 7 Q |y|l+M 2 

It is seen that at the limiting value fi = 0, T(0) is divergent, so that the integral (fTB"]) is divergent as well. The same 
happens for fi = 2, in view of the divergence of T(— 2). However, irrespective of how close to or 2 the label fi > 
is, the integral (|T3"|) is convergent. 

It is interesting to observe that the divergence of the the Fourier integral, as fi approaches or 2, becomes com- 
pensated, if we substitute it back to Eq. (fT2"]) and next consider the limiting behavior of the result: 

2r(l + u)r(-/i)sin^cos^ f°° 
-|A|^ 2 /(x) = — i w V \_ 2 2- / \k\»f(k)e- lkx dk = 

k^fitye-^dk. (14) 



V2^ J- 
Here we use the identity 



r(i + M)r(-/i) = -- — . (is) 

sm ir/i 



In view of the above divergence obstacle, the integral representation (|TU|) is invalid at the boundaries of the stability 
interval. However, the range of validity of its ultimate Fourier version, e.g. the right- hand-side of Eq. (|14[) . can be 
safely extended to the boundary values and 2. 

In fact, it is ——j^ J^°cx> f(k)e~ lkx dk = — d^/dlx^ = (— A)^/ 2 that is commonly interpreted in the literature as 
a definition of the fractional derivative of the fi-th order. On formal grounds this definition encompasses the boundary 
cases (—A) = 1 and —A. 



A. Heavy-tailed case 

Let us first consider a specific pdf (member of the so-called Cauchy family, p(x) = 2/[7r(l + x 2 ) 2 ], such 
that p 1 / 2 {x) = yf j^x 1 an d F° ur i cr transform reads p 1 ^ 2 {k) = e~' fe '. In view of e~ %kx f{\k\)dk = 
2 J °° cos(kx) f(k)dk we have 

/•OO 

V M (x) = -(1 + x 2 ) / fc' 1 e~ fc cos for dk = -(1 + x 2 ) 1 ^ iL T(l + ^)cos[(l + /i)arctana;] , < fi < 2. (16) 
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Expression (fT6|) permits to reproduce easily a number of for different values of p S (0, 2). All V^, p £ (0, 2), derive 
from a pre-defined (quadratic Cauchy) pdf. However, we find most interesting the following three limiting cases. 
Namely, for fi = 1, we get the result derived previously in Ref. [f| 



In the vicinity of the boundary value p = 2 the expression (|16[) yields 



V^ a , 2 (a;) « 2 ^ - L [^(^ - 3 ) arctanx + (3a; 2 - 1) (2 7 - 3 + ln(l + x 2 ))] , (18) 



where 7 w 0.577216 is Euler constant. Hence, for the boundary value /i = 2we get 



, , 2(3x 2 - 1) 

= (19) 



which stays in conformity with a naive expectation that — |A| = A. Indeed, of formal grounds, we readily get (| 19[) 
via Vi{x)p 1 l 2 {x) = +d 2 p 1 /' 2 (x)/dx 2 , as expected. 

For p = 0, in view of arctanx = arccos (1/Vl + z 2 ), from (12) we recover Vq(x) = — 1 which is compatible with 
Vq(x) — — ^!/2 — = — 1 upon an identification |A|° = 1. 



B. Gaussian pdf 

Let us consider the invariant pdf in a Gaussian form p* = — l=e~2^ whose square root p~J 2 has the Fourier image 
(p*) 1/2 (k) = f(k) = ^j^-^= [°° e^e^dx = ^^- e ~ k2 ° 2 . (20) 



Accordingly 



w = -(7)17^)^ / (pnnmre-^dk 



2V^ 



( P *y/2(x)^^ Jo 



Wcosfcx e~ k a 'dk. (21) 



Clearly, at p — 2 we arrive at 



whose equivalent derivation is provided by setting — |A| = +A in Eq. (J7J. Indeed, we verify by inspection that 
v„( x ) - f"W 

The case of p = can be handled as before with the outcomes Vo(x) = —1 and |A|° = 1. To this end we observe 
that (set p, = in Eq. (13)) J °° cosfcx e'^^dk = ^|exp(-^r) while (p*) 1/2 = -^Lexp(-^). 



III. /x-FAMILY OF PDFS FOR A PREDEFINED V(x): GENERALTIES. 

Presently, we shall proceed in reverse. If we choose a priori a concrete potential function V(x), then an ultimate 

1/2 

functional form of an invariant pdf p*{x) (actually pj (x)), if in existence, needs to come out from the compatibility 
condition (|7J|. For any functional form of the confining potential V(x), Eq. ([7]) imposes the following identity to be 

1 /2 

obeyed by an invariant (terminal) pdf (we denote pj (x) = f{x)) 



V(x)/(x) = -|Ar/ 2 /(x), 



(23) 
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where < fi < 2. Remembering that we consider V[x) to be a continuous and bounded from below function (may 
be unbounded from above), we turn over to the standard Fourier transform method with f(k) = -J= J_ ao f(x)e dx 

and fix) = — J= f(k)e~~ lkx dk. Denoting Uk the Fourier image of the right-hand side of Eq. (|2"3"|) . we obtain 



u k = -\k\»f(k). 

Equating Fourier images of both sides of Eq. (|23|). provided they exist, yields 

^ k - k ^f{k')dk'dx. 



If 1 ff 



(24) 



(25) 



In this case, the Fourier image f(k) of a solution f(x) to Eq. ([23]) is defined by following integral equation of a 
convolution type 



/(*) = - 



1 



fc|^V27T 



V{k-k')f(k')dk'. 



(26) 



Here we employ the identity e ™(k-k )^ x _ _ valid in the sense of distributions. 

Our confining assumption implies that the function Vp,{x) typically grows at infinities so that its Fourier image 
may not exist, unless distributionally. The same jeopardy appears in case of a product V(x)f(x), (|2"5)) . Modulo those 
obstacles, we shall demonstrate that, if interpreted in the sense of distributions, the general solution of Eq. ([2"B[) can 
be expressed in terms of (^-function and its derivatives. 

Let us restrict further consideration to even functions V (in this case the terminal pdf is an even function as well) 
and for clarity of presentation assume them to be entire functions. The Taylor series comprise even powers of x only 



vW = v(o) + v'V(o)| + v( 4 »(o)^ + 



4! 



(27) 



The Fourier image of ([27). c.f. ([25]) and ([26]). yields 

1 



V(k) 



'2it 



V(0) + V"(0)-+VW(0)- 
V( W )-^(*) + ^^ W 



c dx 



(28) 



i.e. it has the form of the infinite series of even derivatives of the Dirac S - function. We note that if V(x) is a simple 



(even) polynomial (example: V(x) 



bx ), the above series are finite. 



Accordingly, we end up with the following differential equation of the infinite even order for the Fourier image 
f(-k) = f(k) of pl /2 (x) = f(x): 



v"(o)d 2 /(fc) yW(o)rf 4 /(fc) 

2 dk 2 4! dk 4 

We choose the following initial conditions for Eq. ([2"9")l 



P+V(0) 



/(AO, k>0. 



f(k = 0) = 



f(x)dx = A,f 



(2n-l)/ h _ n.\ _ 



(k = 0) = 0, n= 1,2,3, 



(29) 



(30) 



Note that this imposes an integrability condition on p 1 J 2 {x) on 1Z. Here by /' 2n_1 '(A; = 0) we denote the odd 
derivatives of f(k) at k = 0. 

The integration constant A is not completely arbitrary and should be consistent with the normalization condition 
f 2 { x )dx = p*(x)dx = 1. In view of the Parceval identity we have 

/oo />oo 
f 2 (k)dk = 2 f 2 (k)dk = 1. (31) 
-co JO 

This means that f(k = 0) = A must be compatible with f 2 {k)dk = 1/2. 

One should not expect an easy analytic outcome of the solution of the infinite order differential equation ([2"§]) . In 
most cases of interest the infinite series can be truncated, but generically a numerical assistance is unavoidable. The 
practical strategy of finding (at worst approximately for truncated series and an arbitrary functional shape of V(x)) 
an L 2 (1Z) integrable non- negative ground state of the a priori prescribed semigroup, and thence the terminal pdf 
p*{x), can be summarized as follows. 
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• Expand V(x) in power series. The number of terms in the series should be chosen so as to obtain a sufficiently 
good approximation of the potential. 

• Solve the differential equation (|29[) with initial conditions (|30j). If V(x) is a polynomial function, there are good 
chances to solve this equation analytically. Otherwise we should reiterate to numerics. Check a compatibility 
of f(k = 0) = A with the normalization condition. 

• Analytically or numerically take the inverse Fourier transform to obtain a non-negative function f(x), to be 
interpreted as p 1 J 2 {x). 

• Check whether the obtained f(x), x 2 f(x) and |A|^ 2 /(a;) are absolutely integrable (to guarantee that they 
actually are Fourier transformable). Additionally check that they are square integrable functions and L^R) 
normalize the resultant f(x). 

• Everything is being accomplished under an assumption that | Ap'/ 2 + V{x) is not merely symmetric, but a self- 
adjoint operator in a suitable domain that is dense in L 2 (R). A tacit assumption is that p\' 2 (x) actually is in 
the domain of H. However, this point is quite delicate and involves jeopardies, as we shall see in what follows. 
The operators considered are always symmetric. An issue of their self-adjoint extensions is by no means trivial. 
Only a self-adjoint generator allows to interpret the derived pdf as the terminal (asymptotic invariant) one for 
an associated semigroup dynamics 



IV. CASE STUDY OF THE PDF RECONSTRUCTION FOR V = x 2 /2 AND p e (0, 2) 

We shall consider an exemplary analytic (in part computer-assisted) realization of the previously outlined inverse 
procedure. Our motivations come from Ref. where the Cauchy oscillator has been addressed. 

We recall that by analogy with the familiar harmonic oscillator Hamiltonian (related to the Ornstein-Uhlenbeck 

process) H = —DA + -^j , whose non-negative spectrum starts from 0, we have considered, [Edl, the Cauchy 

oscillator problem in the form JT1/2 = A|V| + x 2 — £0), with £q left unspecified. The compatibility condition (6) 

appears in the form x 2 — £q) p^J 2 = —A |V| p X J 2 . Its Fourier transform reads 

-|A p / + 7H/ = £b/ (32) 
1/2 

where f(p) stands for the Fourier transform of / = pj (x). 

By changing an independent variable p to k = (p — cr)/£, next denoting ip(k) = f(p) with the identifications 
a = £0/7 and £ = (4K/7) 1 / 3 , we may rewrite the above eigenvalue problem (with £q standing for an eigenvalue) in 
the form of the following ordinary differential equation 

^-=2\k\m, (33) 

whose solutions can be represented in terms of Airy functions, c.f. Ref. We note in passing that a slightly different 

scaling was originally employed in Ref. @: C = ( K /27) 1 ^ 3 — > d = \k\ip{k). 

The equation (|33p is in fact a departure point for our subsequent discussion. Quite at variance with the standard 
(text-book) harmonic oscillator intuitions and our own discussion of the Cauchy oscillator spectral problem Q, just 
out of curiosity and as a useful exercise let us consider the compatibility condition for V = .t 2 /2 proper, e.g. without 
any additive counterterms of the form — £q: 

^pl /2 = -\^ 2 p 1 J 2 1 i)<p<2. (34) 

That amounts to assigning the "improper" eigenvalue zero as the bottom eigenvalue of the corresponding spectral 
problem. 

We take Fourier images of both sides of Eq. (f3~4"| to obtain 

1 f°° x \r \ l k XJ 11 d 2 f°° ,. . lkxj ld 2 f(k) 



2W_oo 2 Jy ' 2J2^dk 2 J_ OD Jy ' ~ 2 dk 
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Accordingly, we have 



d 2 f(k) 
dk 2 



2|fc|"/(fc) 



(36) 



whose special [i = 1 case the previous Eq. (|33p actually is. 

To find a solution to Eq. (|36[) . our main idea herefand possibly for an arbitrary potential) stems from the approach 
we have originally adopted for Eqs. (|3"2j) and ([3"3")) . (El Il9j and based on the exploitation of Airy functions. We note in 
passing that both the (Fourier) integral definition of the fractional operator and the solution construction procedure 
involves the Cauchy principal value integral and the continuity (gluing) conditions in the vicinity and ultimately at 
the point £ R. 

We should find the decaying solution of the corresponding differential equation Eq. (I36[) in the fc-spacc, on the 
positive semi-axis (fc > 0), and an oscillatory one on the negative semi-axis (fc < 0). Then we need to shift the 
obtained solution to the right so that the first maximum of the oscillatory part is located at fc = 0. After "chopping" 
the rest of the oscillating part one has to reflect the remaining piece about the vertical axis to get an even "bell-shaped" 
function. 

The obtained fc-space solution should be Fourier-inverted (properly - here we encounter a crucial difference with the 
reasoning of [E[ for the Cauchy oscillator) and squared, while keeping in mind the L 2 (R) normalization issue. This 
procedure is expected to determine the desired invariant (and a candidate for terminal) pdf in the x -space. 

Remark 3: We need to explain how our workings with the +k m shift are to be understood, since (|36[) itself is not 
translationally invariant. Eq. (|36p has the form d 2 f/dk 2 = G(fc)/(fc), where G(k) can in principle be an arbitrary 
function. We execute k — > k + k m everywhere in (f3"o| , so arriving at d 2 / /dk 2 = G(k + k m )f(k + k m ). Obviously, with 
k! = k + k m we have d 2 f/dk' = G(fc')/(fc'), where fc' £ R and the same as previously boundary data at infinities. 
This function is actually Fourier inverted with respect to the fc'-label. To see clearly, why such pre-caution needs to 
be kept in mind, lets us notice that f(x) = / °° f(k')cos(k'x)dk' maps back (|36[) into ([34]) . On the other hand, to map 

([33)1 into the spectral problem for if 1/2 = A|V| + x 2 — Eq), we need to execute the inverse Fourier transform with 
respect to p in ip(k) = ip[k(p)} = f(p). 

To solve ([36]) and next ([34]) . according to the previously outlined procedure, let us consider a pair of equations 



= 2 signfc |fc| M /(fc), instead of the single ([36]) 



d 2 f(k) 



fc > 



= -2(-krf(k), k < 0. 



(37) 



The resultant solutions have different forms for fc > and fc < respectively, [16(. Namely, for fc > we have 



/(fc) = sfk 




C 12 K± 

2q 



'V2 



k" 



;(M + 2), 



(38) 



while for fc < there holds 



/(fc) 



2q \ q 




(39) 



Here J v {x) and N v (x) are Bessel functions and I v {x) and K u (x) are modified Bessel functions, see [131 ■ The 
asymptotics of I„(x) and K v {x) at x —¥ 00 read [l7| 



K v {x) 



(40) 



while as fc — > —00 the asymptotics of the functions J u (x) and N v [x) is oscillatory fl7l |. This means that to obtain a 
localized pdf, we should leave the term with K±_ in (|38p only, so that /(fc) assumes following form 



/(fc) 



C 12 VkK^ (f fc« 



C2 - J h (f |fcT)+^(f 



fc > 
fc < 0. 



(41) 
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TABLE I: Roots k m {n) of Eq. ((45} corresponding to first maximum of oscillatory part of (|42|l for different /i (middle column) 
and normalization constants C(fj,) (right column). 







C(/x) 


0.0 


-0.55536=- -2- 

4v2 


0.597135=^/| (l + f)~ 1/2 


0.2 


-0.621962 


500134 


0.4 


-0.679458 


0.429855 


0.6 


-0.729002 


0.376894 


0.8 


-0.771717 


0.335701 


f.O 


-0.808617 


0.302823 


1.2 


-0.840577 


0.276010 


1.4 


-0.868346 


0.253745 


1.6 


-0.892550 


0.234970 


1.8 


-0.913716 


0.218927 


2.0 


-0.932286 


0.205597 



As the equation (f3"T|) is of second order, we should impose the continuity conditions at k = for a solution and its 
first derivative. We note, that the numerical solution of equation (|3T|) directly involves the value of function and its 
first derivative at k = 0. That implies (see Appendix lAl for more details) 



where C = C12, 



/(*) = c 



1 

2^ 




i^w), fc > 

f [cot ^ J^(u) - 2V„(ii)] , k < 0, 



(42) 



A- 1" 



2^ 



We note here that for /1 = 1 we obtain from 



/(fc) = CVkKi ( '^ki 



3 / 2b 



|i+f 



(2U) 



(43) 



(44) 



known from Ref. (BJ, where the eigenvalue problem for the Cauchy oscillator has been solved, see also Ref. 12 1. 

The next step is to find the position fc m of the first maximum of an oscillating part, next shift the solution to the 
right by k rn , reflect the solution with respect to the y axis and "chop" the rest of oscillating parts. By equating to 
zero the first derivative of an oscillating contribution to (|42[) , after some algebra we get the following equation 



N v -i(u) - cot — J v -i(u) = 0, 



(45) 



where v and u are defined by (|4"3")) . Solutions of this equation can be tabulated, see Table |U The "raw" solutions ([4"2")) 
are shown (along with the position of fc m ) in Fig. [T] 

The normalization condition allows us to fix the admissible values of hitherto unspecified constant C. Namely, we 
have 



C 1 



f{k)dk = 2C 2 / f{k)dk = 2C^ 
'0 



fi(k)dk 



fi(k)dk 



= 1, 



(46) 



where /1 and /2 denote the oscillatory and decaying parts of Eq. (|42j) respectively. This integration can be performed 
numerically and results are reproduced in the right column of the Table |U The final step of the procedure is to 
invert the fc-space solutions to .T-space and square them to obtain the invariant pdf. Except for special /1 cases, this 
procedure can be accomplished only numerically. 
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FIG. 1: Raw solutions of Eq. (|42[) . Figures correspond to /i values. Solution for /i = 1 corresponds to Airy function (|44|) . 
Formal solution for /i = and the position of the first maximum of the oscillatory part k m are shown as an example. 



V. LIMITING (MIS)BEHAVIOR AT THE BOUNDARIES OF (0, 2) 3 p AND THE ZERO (BOTTOM) 

EIGENVALUE ISSUE. 

The stability interval fi £ (0, 2) is an open set. However, since /i can be chosen to be arbitrarily close, respectively 
to or 2, simply out of curiosity it is not useless to address a hitherto unexplored issue, of what is actually going on 
in the limiting behaviors of p I and p, ^ 2, that can be consistently executed on the level of Fourier transforms. 

We note that the operator — |A| M / 2 , as defined by Eq. (fT0|) . is a pseudo-differential (Ricsz) operator and the integral 
there-in needs to be taken as its Cauchy principal value. Hence both in the x and fc-spaces the point G R is 
particularly distinguished. 

On formal grounds, as exemplified in Section II, a naive expectation would amount to literal setting of p = or 
p = 2 instead of the "normal" stability index values p £ (0,2). Then, the operator — |A|^/ 2 would turn over into 
— |A| = —1 and — |A| = A respectively. However, this compelling picture turns out to be problematic, except for 
rather special cases, like those considered before in Section II. 

We recall that the core of our solution method lies in Fourier transforming the compatibility condition (JJJ 

^ P V 2 = -\A^VJ 2 ^^ = 2\krf(k) (47) 

and actually solving the obtained differential equation in the fc-space. 

Given the solution, the unsettled problem is the validity of the inverse Fourier transform of the outcome. That 

means that not only f(k) itself, but also \k\^f(k) and specifically d need to be Fourier invertible to reproduce 

1 /2 

the cc-space version of the problem. One should as well be aware that with a resolved pj , we must have secured the 
validity of the Fourier transform, that is explicit in Eq. (|47l) . 
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A. Case of /i = 



We seriously address a Fourier transformed equation (j47|) and employ again Eqs. (|36|) and (|37|) . while blindly setting 
/i = there- in. The "raw" solution of (p?T|) can be either obtained from (|4"Tj) or explicitly from (|57|) . It has the form 

f C 12 e" fcV 2, fc > 

/(fc) = ^ (48) 
[ C 2 i cos kV2 + C22 sin fc\/2, fc < 0. 

The continuity condition at fc = for /(fc) reads 

C 12 = C21 (49) 

and for the derivative 

C22 = — C12 = — C21, (50) 

so that (we set C12 = C) 

{e" fcV2 , fc > 

(51) 
cos k\[2 - sin fc^2 = \/2 cos (fcV2 + f ) , fc < 0. 

The first maximum of oscillatory part is located at 

fc m = ^= ps -0.555360367, (52) 

4V2 

in accordance with Table |U The "raw" solution (jSTj) is shown (along with the position of k m ) in Fig. [T] 
Now we shift the whole solution to the right and " chop" the unnecessary piece of an oscillatory part 

V2cosfc\/2, < fc < -k m 

(53) 

e -(fe+fc m )V2 k > _ k 



The normalization condition (|46[) reads 



2C* 2 



cos 2 fcV2<ifc + / e~ 2 ^ k+km) dk 



(54) 



so that C = 4 , 7 ps 0.629325. This normalization coefficient is different from that in Table U because the 

vVHt 



transition from Bessel functions of index 1/2 to elementary ones introduces an auxiliary coefficient y / 7t/2 3 / 4 . Thus, 
we have C t 2 able =(2v / 2/7r)C 2 = (2/7r)/ (l+7r/4). This minor difference is insignificant with respect to the normalizability 
of the pertinent eigenfunction, as an overall coefficient before /(fc) is just C, §54 
Now we invert the Fourier transform to get 



2 f°° 2 
fix) = == / f(fc) cos kx dk = C\ — 

K ' Jo V 7T 



V2 / cos fc^2 cos fcx dfc + e fe " iv/2 / e fcV2 cos fcx dfc 



= 5 r ^ sm — 7= — v2cos — = G = 



7T x 4 - 4 V 4\/2 4\/2 / ' y/2Jl 



(55) 

We note that at the point x = ±\/2 function /(x) looks divergent. However, since the terms in the numerator yield 
the same zero as in the denominator, the potentially dangerous divergence is removed when one approaches x = ±y/2. 
The function f(x) decays at spatial infinities as 1/x 3 and is positive for |x| < 5. For |x| > 5 we encounter oscillations 
so that both zeroes and negative values are developed. Thus, it is clear the function f(x) cannot be interpreted as 
an (arithmetic) square root of a probability density p*(x). By this reason alone the obtained solution is incongruent 
with the x-space version of the compatibility condition. 

More than that, an inversion of the Fourier transformed equation Eq. (|4"T|) back to the x-space leads to a con- 
tradiction, since ^-/(x) 7^ —f(x). The roots of this discrepancy are obvious. The function (|55|) at spatial infinities 
behaves as f(x) ~ 4, sin j^=. Consequently, x 2 f(x) ~ ^ sin is not absolutely integrable on R and thus not Fourier 
transformable. 
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B. 



Case of /i = 2 



E )p 



1/2 



1, LO 2 



= 0, where Eq = 1/2 is the lowest 
where m is a particle mass and uj is 



Concerning /i = 2, we recall the eigenvalue equation (—A + ■ 
eigenvalue for a quantum harmonic oscillator in the units h 2 /2m 
an oscillator frequency. That gives rise to the Gaussian ground state function as a celebrated textbook result. 

From this point of view the adopted (Section IV) form of the compatibility condition (possibly, with a notable 
exception of (i = 1) might seem puzzling. Specifically, the standard eigenvalue equation for the harmonic oscillator is 



1/2 



demand, expected (per force) to correspond to the case of fx = 2 in 



plainly incompatible with the (—A + \)p* 
the considerations of Section IV, c.f. Eq. l47l 

Nonetheless, the latter equation seems to admit a non-Gaussian heavy-tailed solution, 
(|38|) -(|42 |) . In particular, an asymptotic (k — > oo) behavior of 



which is derivable from 



K 1/2q (— kA ~k-i' 2 exp 



(56) 



with q = (fj, + 2)/2, /i G (0,2] yields an asymptotics in the fc-space for (j, = 2, cf. (|42|) . as follows: f(k) ~ 
(1/Vk) cxp[-fc 2 /V2]. That secures an absolute convergence of the inverse Fourier integral. 

We have no sufficient analytic tools in hands to analyze various features of the ultimate x-space solution. However, 
the performed numerical (computer assisted) analysis proves that both f(x) and x 2 f(x) are absolutely intcgrable. 
The same occurs for [x 2 f(x)] 2 . The pertinent f(x) = p X J 2 (x) has an inverse polynomial decay at infinities, decaying 
faster than 1/x 3 . In the vicinity of x — the obtained function mimics a Gaussian. Those properties are visualized 
in Figs. 2 and 3. 

The annoying point that should be clearly spelled out in connection with ou derivation procedure may pertain to 
the point G R and the fact that we have actually (albeit tacitly) defined the Laplacian not on the whole of R, but 
on {i?\0}. A subsequent comment plj indicates where the peculiarity of our solution may possibly be rooted. 

Comment 1: To make a clear distinction with the text-book quantum mechanical reasoning for the harmonic 
oscillator problem, let us point out that we have actually obtained appears to be related with the (normally discarded) 
"improper" eigenvalue n = —1/2 in the familiar E n = ftw(n + 1/2) formula. Here we emphasize the text-book 
restriction upon the harmonic oscillator eigenvalues n > —1/2, whose interpretation is: "all states corresponding to 
n < —1/2 need to vanish identically", |2(j. A sufficient condition to this end is that the lowest energy eigenfunction 
vanishes under the action of the annihilation operator. This actually selects n = as the bottom spectral label and 
n G (AT (JO). Our reasoning demonstrates that it is not a necessary condition, and actually, an "improper" eigenvalue 
n = —1/2 is admissible as a legitimate spectral label in the harmonic oscillator problem. However, in this case a 
specific choice of the boundary data (implicit in the construction procedure but hitherto nor explicitly spelled out) 
appears to enter the game. That may have an effect on the spectral features (including that of the proper domain 
definition and the self-adjointness issue) of the original problem, see e.g. the next Comment. 

Comment 2, [2l|: Generically, the Laplacian — A, if defined on a Schwartz space S(R) (fast decaying functions, 
continuous and with continuous derivatives), is non- negative and essentially self-adjoint operator . Therefore its 
closure, which is a self-adjoint operator, is non-negative as well. Its continuous spectrum extends from to +oo. 

However if we define the Laplacian on 5({i?\0}), the operator remains non-negative but ceases to be essentially self- 
adjoint. Its deficiency indices read (2,2), implying an existence of a two-parameter family of self-adjoint extensions. 
There is a well known in the literature example of such extensions determined by the boundary data /' (+0) — /' (— 0) = 
a/(0). On the operator level so point-wise restricted operator — A a may be symbolically represented —A + aS(x), 
where 6(x) stands for Dirac delta distribution. 

For a = we have a "normal" Laplacian. The case of a — >• +oo sets Dirichlet boundary condition at (an 
impenetrable barrier). If a > 0, the operator —A a has a continuous non-negative absolutely continuous spectrum. 

The case of a < is particularly interesting, since then — A a has one non-degenerate isolated negative eigenvalue 
— a 2 /2 with an eigenfunction ( | a| x / 2 exp(— |ax|/2), see e.g. [UJ. The remaining part of the spectrum begins at and 
is absolutely continuous up to +oo. 



An operator —A 



defined on 5({i?\0}), is regarded as positive and symmetric, but is not essentially self- 



adjoint. The deficiency indices are likewise (2,2). Assuming the previous boundary conditions we turn over to the 
one-parameter family of self-adjoint extensions — A a + x 2 , assuming likewise a < 0. 

The choice of a = leads to the standard harmonic oscillator problem, while for negative a there appears a fairly 
non-standard spectral solution described in detail in Rcf. 



22j. The problem has exactly one negative eigenvalue plus 



a positive discrete spectrum, of the form e = (2v + l)/2. 
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-4-2 2 

X 



FIG. 2: f(x) = pl^ 2 (x) for /i = 0, 0.3, 1, 2. The logarithmic scale is employed to better visualize the behavior on the "tails" 

Given negative a, the eigenvalues e come out from a transcendental equation v — oF(l — i//2)/r(l/2 — i//2) = 0. 
The value e = appears for v = —1/2 i.e. for a = — r(3/4)/T(5/4). The corresponding eigcnfunction shows up a 
~ x -1 / 2 exp(— x 2 /2) decay at +oo, |22| . 

VI. CONCLUSIONS 

To conclude, here we have presented an outline of the general formalism for how to find an invariant pdf for a 
pre-defined semigroup potential and all stability index values < \i < 2. We have considered the generic case of 
symmetric even pdfs, for which a potential V(x) is an even function. That was dictated by known spectral properties 
of fractional operators associated with a symmetric stable noise. We have reduced the problem of finding a terminal 
pdf to that of finding a solution of the ordinary differential equation with infinite number of terms in momentum 
space. Such differential equation, even if hard to handle analytically, can rather easily be solved with a numerical 
assistance. For polynomial potentials the number of terms becomes finite and the pertinent equation can be solved 
analytically. 

The outlined procedure has been explicitly performed for the family of Levy stable oscillators, with a common 
quadratic semigroup potential V(x) — x 2 /2. In this case, a solution of the corresponding differential equation in 
the fc-space has been obtained by employing a suitable continuity procedure at k = 0. After Fourier-inversion and 
squaring the result, this yields an ultimate functional form of the desired invariant (and possibly terminal) pdf of the 
jump-type process, for arbitrary pi <G (0,2). 

We have analyzed a limiting behavior of solutions in the vicinity and at the boundaries [i = and 2 of an open 
stability interval (0,2). In the case of /i = we have derived an explicit analytical form of the solution in k and (by 
Fourier inversion) in the x-spaces. The pertinent function shows a nontrivial oscillating behavior and thus is definitely 
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not a square root of any pdf. However, with an explicit analytic form of f(x) one easily proves that the x 2 f(x) is 
not-Fourier transformable, hence the crucial connection between the x-space and k-space solutions, Eq. (|T7|) has been 
lost. 

1 li 

We have shown that for fj, = 2 a positive pj (x) is obtained. This invariant solution definitely corresponds to 
the "improper eigenvalue" e = of the energy operator. To comply with the common knowledge about the spectral 
properties of the operator — A + x 2 , we find that the (hitherto neglected) presence of the special boundary data has 
been enforced. Then only the self-adjoint extensions of the corresponding non-negative operator may be introduced, 

1 /2 

and pj (x) may receive a consistent interpretation as an asymptotic (terminal) target in the semigroup evolution. 

This particular case in turn, indicates the main jeopardy hidden in the presented reconstruction problem. Once an 
invariant pdf has been settled, it is necessary to identify a consistent self-adjointness domain for the involved operator 
H . Then only an invariant pdf may be interpreted as an asymptotic (terminal) one for a well defined semigroup 
exp(-tH). 

A problem of reconstructing a self-adjoint operator from the knowledge of its ground state is not trivial at all, see 
e.g. for an investigation of [23| in the diffusion-type framework. In the harmonic oscillator problem of Comment 
2 it is known that a corresponding self-adjoint Hamiltonian shares odd-parity (label n = 1,3,5,...) part of the 
spectrum with a "normal" harmonic oscillator, while an even-labeled part is substantially different (that pertains to 
the associated eigenfunctions as well). 

Acknowledgement: We are willing words of gratitude to Professor Witold Karwowski for enlightening discussion 
and explanation of the role of boundary data in differentiating between the " normal" and " abnormal" (zero bottom 
eigenvalue) harmonic oscillator spectral problems. 
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Appendix A: Continuity conditions at k — 
From (|4"Tj) . the functions at k = read 

k=0- 

The derivatives at k = 

C 2 i[\/*J„(tO]fe=o + C 22 [\fkN v {u)]' k=Q = CiatVAJ^^Uo. 



(Al) 



(A2) 



Such forms of (|A1[) and (| A2|) are dictated by the following asymptotic expansions of Bessel functions near k = in 
variables (|4"5j) 



X-„(u) 



JV„(u) 



J„(u) 



r(-u) 



Vk 



V2 
H + 2 



COS 7W r( — f) 



V2 



M + 2 



V2 

r>) 

•K\fk 



M + 2 



r(i + i/) 



M + 2 



(A3) 



Eq. (| A3|) means that {\fkJ v {u)]k=a = 0. For reference purposes, the derivatives like WkN v {u)]' k=Q (we first multiply 
by y/k and then differentiate) read 



[VkN v (u)}' k=0 = - 



cos wis T(—v) 



V2 
Ai + 2 



, [VkK„(u)}' k=0 = 



r(-i/) 



M + 2 



V2 



At + 2 



Substitution of values of functions and derivatives into (|A1[) and (|A2|) yields 



7T 7T 7W 

C22 = — ■^■C'12, C21 = 2^ 12 co * "2"' 



(A4) 



(A5) 



which, after employing the identity T(l + /i)T(— fj,) = — 7r/ sin(7r/x) gives rise to Eq. (|42|). 



[1] R. Metzler and J. Klafter, Physics Reports, 339 1, (2000) 

[2] I. Eliazar and J. Klafter, J. Stat. Phys. Ill, 739, (2003) 

[3] P. Garbaczewski and V. A. Stephanovich, Phys. Rev. E 80, 031113, (2009) 

[4] P.Garbaczewski, V. Stephanovich, Phys. Rev. E 64, 011142, (2011) 

[5] P.Garbaczewski, V. Stephanovich, Physica A, 389 (2010), 4419 

[6] P. Garbaczewski, V. A. Stephanovich and D. Kedzierski, Physica A 390, 990, (2011) 
[7] D. Brockmann, I. M. Sokolov, Chemical Physics, 284, 409, (2002) 

[8] E. B. Davies, Heat kernels and spectral theory, Cambridge Univ. Press, Cambridge, 1989 
[9] J. Bertoin, Levy processes, Cambridge Univ. Press, Cambridge 1996 
[10] K. Kaleta and T. Kulczycki, Potential Analysis, 33, 313, (2010) 

[11] K. Kaleta and J. Lorinczi, Analytic properties of fractional Schrddinger semigroups and Gibbs measures for symmetric 
stable processes, larXiv: 1011.27131 (2010) 

[12] J. Lorinczi and J. Malecki, Spectral properties of the Cauchy massless relativistic harmonic oscillator, arXiv: 1006.3665 
(2011) 

[13] R.Carmona, Commun. Math. Phys. 62, 65, (1978) 

[14] R. Carmona, W. C. Masters and B. Simon, J. Funct. Anal. 91, 117, (1990) 

[15] B. Dybiec, I.M. Sokolov, A.V. Chechkin Journal of Statistical Mechanics: Theory and Experiment (2010), p07008. 

[16] A.D. Polyanin, V.F. Zaitsev Handbook of Exact Solutions of Ordinary Differential Equations. CRC Press, 1995. N 2.1.2.7. 

[17] Handbook of Mathematical functions. Ed. M. Abramowitz and I. Stegun. National Bureau of Standards, 1972. 



15 



[18] A. A. Dubkov, B. Spagnolo, V.V. Uchaikin, Internat. J. Bifur. Chaos 18 (2008) 2549. 
[19] R.W. Robinett Amer. Journ. Phys., 63 (1995), 823. 

[20] R. L. Liboff, Introductory Quantum Mechanics, Holden-Day, San Francisco, 1980 
[21] W. Karwowski, private communication 

[22] J. Viana-Gomes and N. M. R. Peres, Eur. J. Phys. 32, 1377, (2011) 
[23] R. Vilela Mendes, J. Math. Phys. bf 27, 178, (1986) 



